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!>■ ■ Abstract 

o 

We study the QCD scattering amplitudes for qq — ► gg and — > ggg where q is a massive 
fermion. Using a particular choice of massive fermion spinor we are able to derive very 
compact expressions for the partial spin amplitudes for the 2 — > 2 process. We then 
investigate the corresponding 2^3 amplitudes using the BCFW recursion technique. For 
the helicity conserving partial amplitudes we again derive very compact expressions, but 
were unable to treat the helicity-fiip amplitudes recursively, except for the case where all 
the gluon helicities are the same. We therefore evaluate the remaining partial amplitudes 
using standard Feynman diagram techniques. 



1 Introduction 



In the last two years there has been dramatic progress in the calculation of multi-particle 
scattering amplitudes in quantum field theory. Following Witten's pQ discovery of a connection 
between QCD amplitudes and twistor string theory, a calculational technique |2] was found 
which has come to be known as 'the CSW construction'. It amounts to an effective scalar 
perturbation theory, in which MHV amplitudes are elevated to the status of vertices, connected 
by scalar propagators. This scheme found wide application [3J HJ El EE El IHj , though it turned 
out that there was an even more efficient way to calculate scattering amplitudes. Britto, 
Cachazo, Feng and Witten found a recursion relation which, by shifting momenta, takes 

advantage of the analytic properites of tree amplitudes. Use of the BCFW recursion relation 
led easily to very compact expressions. Originally applied to purely gluonic tree amplitudes, 
the recursion has since been extended to include fermions ^2 gravitons ^3] and loop 
amplitudes ^3 El QUI • As well as perhaps giving hints of as yet unknown mathematical 
structure beyond the Standard Model, these developments are potentially important for the 
calculation of Standard Model backgrounds at colliders such as the LHC. The more accurately 
the relevant cross sections are known, the higher the discovery potential of the machine will be. 

One area of very recent progress is the calculation of amplitudes involving massive fermions. 
It was shown in Ref. ^H] how to generalize Supersymmetric Ward Identities ^H] to include 
massive particles. In this way, different amplitudes involving fields belonging to the same 
supersymmetric multiplet are related by a rotation. For instance [2*U1 . amplitudes involving 
quarks and gluons are related by SWIs to amplitudes involving scalars and gluons, and these 
have been calculated in Ref. [2^. The off-shell Berends-Giele j22j recursion has also proved 
useful (221- Tree amplitudes with massive fermions are required as input within the unitarity 
[21| method to calculate one-loop amplitudes, and to this end Ref. provides four- and 
five-point amplitudes with D-dimensional fermions, calculated using BCFW recursion. 

The recursion relations were extended in Ref. [211 to include massive fermions, and in [2Z] 
four-point amplitudes involving two massive quarks and two gluons were calculated. The topic 
of the present work is to extend this to five-point amplitudes, using a treatment of massive 
fermion spinors introduced some twenty years ago in Ref. j2Hl- We first outline the particular 
spinor helicity method we will use, and then we show that 2 — > 2 scattering processes in QCD 
can be written in a form which is ideally suited for use in BCFW recursion. In Section 4 we use 
the recursion relations to derive compact expressions for certain qq — > ggg partial amplitudes. 
Those partial amplitudes which we could not treat recursively are evaluated using Feynman 
diagrams in Section 5. Finally, we present our conclusions. 
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2 Spinor Products 



For many years amplitudes involving massless momenta pi and pj have been expressed in terms 
of spinor products, 

[ij] = u + (pi)u~(pj) and (ij) = u~ (pi)u + (p^ . (2.1) 

In this way amplitudes find their simplest expression. The spinors in question can be thought 
of either as two-component Weyl or 4-component Dirac spinors. Numerical evaluation of such 
amplitudes involves the use of the standard formulae for the spinor products in terms of the 
momentum 4- vectors. Following let us first introduce two 4- vectors k and k\ such that 



k ■ k = , fci • fa = -1 , k ■ k x = 0. (2.2) 
We now define a basic spinor U-(k ) via 

u-(k )u-(k ) = ^-to, (2.3) 
and choose the corresponding positive helicity state to be 

u + (k ) = t lU -(ko). (2.4) 
Using these definitions it is possible to construct spinors for any null momentum p as follows: 

u»W = (2.5) 
V2p ■ fc 

with A = ±. Note that this satisfies the massless Dirac equation fu{p) = 0, as required. We 
can now simply evaluate spinor products. For example, 

[ij] = u + (p l )u-(p j ) (2.6) 
u-(fco) piljfj u+(k ) 

(2.8) 
(2.9) 



a/4 (pi ■ fc )(pj • k ) 
Tr(^ % $ # Vi; 



a/4 ^ • ko)(p j ■ k Q ) 

(Pi ■ k Q )(pj ■ hi) - (pj ■ k )(pj ■ fej - ie^pyk^pPkl 
y / (p l ■ k )(Pj ■ k ) 



A similar expression is obtained for the angle product (ij). The arbitrary ko and k\ can now 
be chosen so as to yield as simple an expression for the product [ij] and (ij) as possible, to 
facilitate numerical evaluation of the ampitudes. The choice 1 

fc = (1,0,0, 1) (2.10) 
^ = (0,0,1,0) (2.11) 



lr The notation is = (fc°,k). 
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is a good one, giving the familiar 



ifi + iPi) 



p] ~ Pj 
.Pi - Pi 



(P v i + iPi) 



Pi - Pi 



P° 



Pi 



(2.12) 



2.1 Massive spinors 

To evaluate spinor products involving massive spinors, we need to find a definition analogous 
to ()2.5|) . One possibility is that outlined in 



u\p) 



(f + m)u x (k ) 
V 2 P ■ k 



(2.13) 



which satisfies the massive Dirac equation, (fl — m)u x (p) = 0. The m in (J2.13j) is positive or 
negative when u x (p) describes a particle or antiparticle respectively. This definition has the 
virtue 2 of being smooth in the limit m — > 0. We will use ()2.13|) to evaluate products involving 
massive spinors. 

It is easily seen that the familiar [. .] and (. .) products take the same form for massive spinors 
as they do for massless ones. Explicit mass terms drop out due to various trace theorems. 
However, the product of like-helicity spinors is now non-zero: 



(V) 



M ± (p i )u ± (p j ), 

j 

Pj ■ k 



m 



rrii 



Pi ■ k 

z 
Pj ~ Pj 

P\ ~ P\ 



+ i<r+ j 



J: 



(2.14) 

(2.15) 
(2.16) 



where in the last line we have used k as given in (|2.1L)|) . Note that the like-helicity product is 
the same whatever the helicity of the spinors involved, and that we use a round bracket as a 
shorthand notation for it. 

We have been using the word 'helicity' to refer to the spin projection of massive fermions, 
but in fact this is only justified if the projection is onto the direction of the momentum vector, 
and it is not obvious that this is the case. There exists a unique polarization vector, though it 
depends on the arbitrary /c , 

1 ^ m 2 . „\ . 

(2.17) 



1 

m 



P ■ k 



The spinors (|2.13J) satisfy 



;i - xr 



0. 



(2.18) 



We see that besides the momentum p there is an additional contribution to the polarization 
vector proportional to k . Suppose we have an anti-fermion i and fermion j in the initial state 



2 Care is needed when p ■ ko also vanishes in this limit, as we will discuss later. 
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and they approach along the z axis, in the positive and negative directions respectively. If we 
choose k to be a unit vector in the z direction, i.e. 



k = (1,0,0,1), 



(2.19) 



then for momenta 3 



Vl = (E,O,O,0E), 
Pi = (E,0,0,-PE), 



(2.20) 
(2.21) 



we have the following polarization vectors: 






-(-E(3,0,0,-E), 
rii 

-(EP,0,0,-E). 



(2.22) 



(2.23) 



If we recall that m« is negative because i is an antiparticle, then we see that each polarization 
vector points in the same direction as the corresponding momentum, so that the spinors u x (p) 
are indeed helicity eigenstates for this choice of k^. However, choosing k to be parallel to one 
of the particle's momenta results, in the massless limit, in the denominators of products such as 
that in 1)2.16)1 vanishing. By being careful to take the limit algebraically this does not present a 
problem. 4 But it should be noted that in such cases products like (ij) do not necessarily vanish 
in the massless limit. We can sidestep this issue by choosing a different k , though we could 
not then talk of the helicity of the fermion. 

2.2 Example: qq —> gg 

To demonstrate the use of the massive spinor products described in the previous section we 
calculate the helicity amplitudes M XlX2XiXi for the simple QCD process q Xl (p\) q X2 (p2) — ► 
g X3 {p3) 9 X4 (pa)- The Ai,A2 = ± labels on the quarks refer to their spin polarizations in the 
sense already indicated. If we choose k^ appropriately then they can be thought of as helicity 
labels. We will evaluate the partial (colour) amplitudes for the above scattering process, i.e. 
we consider contributions only from those diagrams with a particular ordering of the external 
gluons. The full colour-summed amplitudes can then be recovered by inserting appropriate 
colour factors, as described in Appendix B. 

We first consider the M + ~ H partial amplitude, for which there are two Feynman diagrams, 
shown in Figure 1. We will express them in terms of massive spinor products. For the slashed 



4 If we take = (1, 0, sin#, cos#), then for the momenta defined in i|2.20|) and (|2.21() . with rrij = —mi = m, 
we have (ij) — — 2m/3cos#(l — 1 cos 2 6*) -1 / 2 . Thus (ij) ~ 0(m) as m — > except if 9 = 0° when (ij) ~ O(E). 
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(a) Diagram A (b) Diagram B 

Figure 1: Diagrams contributing to the colour-ordered partial amplitude for the process 

gluon polarization vectors we use 

ffak) = ^ "4^ + (p)+ M -(p)^-W > (224) 

r(p,ib) = ^ M + (P)"4fc)+^W"-(P) , (2.25) 

where A; is a (null) reference vector which may be chosen separately for each gluon. Different 
choices of reference vector amount to working in different gauges. The choice k 3 = p 4 and 
ki = ps is particularly convenient in this context, as Diagram B vanishes in this gauge. We 
have for the other diagram 

iPl) ~W {P2-PZ?-^ ~W ^ ( } 

which simplifies easily to 

_+/ +, , u + ( Pl )[u-(p 3 )u-(p±) +u + (p 4 )u + (p 3 )]u-(p 2 ) 

« R J^r P4 , 2.27 

4 p 3 ■ Pa Pi ■ Pi 

so that 

M + - + - = [3 |2|4>MM±MM. ( 2 . 28) 
4 P3 ■ Pi Pi ■ Pi 

As promised, we are left with an expression for the amplitude in terms of vector products and 
massive spinor products. 

We next consider the other M AlA2+ ~ amplitudes. It is interesting to note that these are 
directly obtained from the M + ~ + ~ amplitude simply by changing the type of certain brackets. 
Thus 

M+++ - = [131(42) + (14)(32) 

4p 3 -p 4 P4-Pi 

M-+- = |3 |2|4)MM±MS, (2 .30) 

4 P3 ■ Pi Pi ■ Pi 

M-++- = [3|2|4)MM±MM. ( 2 . 31 ) 

4 Ps ■ Pi Pi ■ Pi 
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Those amplitudes where the gluons have helicities (— +) can be obtained directly from the 
ones above by complex conjugation. 

Let us now examine the case where the gluons have the same helicity. By direct calculation 
we find 

M---^[43]"^MM. (2 .32) 
L J (34) 2 2 p 4 - Pl 1 ; 



from which we deduce 



M++++ = (2.33) 

M + ~ ++ = 0, (2.34) 

r i<13>(4 

m[34](12) 



M _ +++ = (13)W - (14X32) 3 

(34) 2 2 p A ■ pi 



(34) 2 Pi - Pl - 



(2.36) 



where in the last line we have used the Schouten identity. The amplitudes with two negative 
helicity gluons are obtained via complex conjugation. There are several interesting things to 
note about these results. First, the amplitude M + ~ ++ vanishes (for any choice of k ) because 
of the identity 5 (13) (42) — (14) (32) = 0. Second, when fc is parallel to the line of approach of 
the fermions (i.e. when we work in the helicity basis) then the product (12), and hence M h++ , 
vanishes. 

We have verified that when squared and summed over spins and colours, the set of 2 — > 2 
scattering amplitudes given above matches the well-known result (see for example Ref. [29 ) 
calculated using Feynman diagrams and 'trace technology', namely 

££M 2 = ^(^-g(n 2 + x! + ,-4), (2.37) 

colours spins 

where 

2p x ■ p 3 2pi • p 4 4m 2 2 
ti = , r 2 = , p = , s = (p 1 +p 2 ). (2.38) 

s s s 

Finally, the m — > behaviour of the spin amplitudes can easily be read off from the expres- 
sions given above. For example, if E denotes the typical scale of the 2^2 scattering 6 , then in 
the m/E — > limit we have 





-+±=F 


M~ 


-±T 


~ 0(1), 


M H 


— ±=F 


M~" 


fT± 


~ 0(m/E), 


M A 


-+±± 


M- 


-TT 


~ 0(m 2 /E 2 






M~ 


f++ 


~ 0(m/E), 


M H 


— ++ 


M~ 


\ 


= 0. 



(2.39) 



5 See Appendix A for a list of identities and notation. 
6 We explicitly exclude zero angle scattering. 
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Note that in deriving these results we have assumed that k is not directed along any of the 
particle momenta, so that all (ij) spinor products are 0(m) in the m — > limit, and (ij), [ij] 
products are 0(E). If on the other hand we choose the (fermion) helicity basis by taking k Q in 
the direction of (say) pi, then (|2.39j) becomes 





M~ 


+T± 


~ 0(1), 


M+- ±± , 


Af - 


+TT 


= o, 


M+ +±T , 


M~ 


-=F± 


~ 0(m/E), 


M++++, 


M~ 




~ 0(m/E), 


M—++, 


M + 


H 


~ 0(m 3 /£ 3 ) 



3 BCFW Recursion 

In Ref. jU] Britto, Cachazo and Feng introduced new recursion relations for amplitudes involving 
gluons. The recursion involved on-shell amplitudes with momenta shifted by a complex amount. 
Later [TOJ , the same authors with Witten gave an impressively simple and general proof of the 
recursion relations. They have since been successfully applied to amplitudes involving fermions 
[TH IT2*] and gravitons ^3] • Risager [3U| has demonstrated how they are related to the earlier 
'MHV rules', providing a proof of the latter simply using BCF recursion. There has also been 
much progress at 1-loop level |HJ EES UH CZj, which has dovetailed nicely with the earlier 
unitarity work [24J. 

We begin by choosing two (massless) particles i and j whose slashed 7 momenta we shift as 
follows, 

where 

ifl = u + (pj)u + (pi) + u~(pi)u~(pj) (3.2) 
is such that both pi and pj remain on-shell. Using the familiar spin-sum condition, 

J2 U \P) u\p) = f (3.3) 

A 

we can re-express the shift (j3.1j) as a shift of spinors: 

+ (3.4) 
u~ (pi) + z u~ (p^ (3.5) 
u + (pj) - z u + (pi) (3.6) 
u~(pj) - z u'{pi). (3.7) 

7 f = l^Pn 



U + (Pi) -> u + (pi) = 

u~{pi) -> u~(pi) = 

U + (pj) -> w + (p}) = 

u"(pj) -> w~(p5) = 
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In the Weyl spinor notation we are shifting Aj and Xj. For massless particles, Dirac 4-spinors 
are effectively two copies of a Weyl 2-spinor, hence the four shifts of (|3.4j) - (j3.7|) . Notice that 
there is no symmetry between i and j — they are treated differently. 

The amplitude is now a complex function of the parameter z. What the authors of [TO] 
showed was that we can use the analytic properties of the amplitude as a function of z to glean 
information about the physical case z = 0. In particular, we get a recursion relation, which can 
be stated as 

A -= E J2 A ^,P- S ) M-P 3 ^)- (3-8) 

partitions s ^ 

where the hatted quantities are the shifted momenta. In fact, this is only valid if the helicities 
of the marked particles are chosen appropriately. The crucial property which must be retained 
if f!3.8|) is to hold is that the shifted amplitude must vanish in the limit z —>■ oo. There are 
rules |1U[ ITU 123 E7J detailing which marking prescriptions are permitted in different cases. For 
our purposes, we will be on safe ground if the shifted gluons have helicites (h\ h?) = (+, — ) or 
(±,±)- 

This method of calculation is particularly efficient because much of the computational com- 
plexity encountered in a Feynman diagram calculation is avoided since the lower point ampli- 
tudes Al and An can be maximally simplified before being inserted in ()3.8j) . 

The sum is over all partitions of the particles into a 'left' group and a 'right' group, subject 
to the requirement that particles i and j are on opposite sides of the divide. The sum over s 
is a sum over the spins of the internal particle. Each diagram is associated with a particular 
value for the complex parameter z, which can be found via the condition that the internal 
momentum P is on-shell. Note that P is always a function of z because of the restriction that 
the marked particles i and j appear on opposite sides of the divide. 

One useful point to note in practice is that three-point gluon vertices vanish for certain 
marking choices. In particular, for the j side of the diagram a gluon vertex with helicites 

(+ H — ) vanishes, as does the combination ( h) on the i side. This was pointed out in 

Ref. 0. 

We will be concerned in this work with the process qq — *> ggg, and so will encounter recursive 
diagrams connected by an internal fermion, the propagator of which is, in this formalism, the 
same as that of a scalar. Following Ref. [27], we 'strip' fermions from the lower point amplitudes 
which feed the recursion and write 



^ = E J2MP,n ^r^r 1 M-P*,Pi), 0.9) 



"p 

partitions s 



E MP,n A R {-P\v 3 ). (3.10) 



rrvp 

partitions 



where P* shows that the amplitude has been stripped of this external spinor wave-function. 
By way of example, let us reconsider the process qfq^ — > otd^- We mark the gluons such that 
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2 



1 



1+ 



5+ 



(a) Diagram A (b) Diagram B 

Figure 2: Recursive diagrams contributing to q + {pi)q~{j>2) — ► g + (P3)g~(p<i)g + (j > 5)- 

i = 3 and j = 4. Then there is one recursive diagram, 

w (Pi) — 7^- 7 ^ j — (P2)- (3.11) 

V2 VP2 - P3> - m v 2 

With the shifts we have chosen, the hats on the polarization vectors can be removed. The 
shifted part of the internal propagator is killed by either of the polarization vectors. So in 
fact all the hats can be removed in (j3.ll)) . which is then identical to the Feynman diagram 
expression ()2.26|) . 



4 qq —> 3g from BCFW Recursion 

The four-point amplitudes we derived in Section 2 are in such a form that it is trivial to strip 
a fermion off in the manner described above. This means that they are particularly convenient 
for use in BCFW recursion. Consider the process qfq^ — > gtg^gt, for which there are three 
recursive diagrams, shown in Fig. 2. We choose the marking prescription % = 3, j = 4. 

The two diagrams with internal gluons both vanish, due to the vanishing of M + h+ and 
the vanishing of the (+ H — ) gluon vertex with the shifts we have chosen. For the remaining 
diagram we use 

M + + = _ [4 | 1 | 3) (1»)M + [14](32) (41) 

4P3'P4P4- Pi 

and strip the fermion u~(p2), leaving 

M+-+ = -[4|1|3> u+( Pl ) v + (P^ + (P*)+v-(P^-(Pz). (4 . 2) 

4 p 3 -Pa Pa -Pi 

After the appropriate relabelling this can be used in Diagram A, which can then be written 

A = -[5|1|4> fi+fo) , ( ^-\ +m \ 

^Pb-PiPb- Pi [P2 - Ps) 2 - m 

x ^-u-{p 2 ). (4.3) 
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Due to our choice of marking, all the hats in the numerator can be removed. The shifted part 
of the propagator is killed by the gluon polarization vector. We are left with 



[5|1|4> 



m(14)(42)[53] + [15] (42)[3|2|4> - (14) [32] [5|1|4> 



8p 5 'PlP2- P3 Pi ■ P5 (43) 

We can work out z from the requirement that P 2 = (p2 — P3) 2 = m 2 , and find 

-2 p 2 ■ P3 



The product p 4 ■ p 5 is then 



Pa -P5 



[3|2|4> 



(P4 - z r])-p 5 

, P2-P3 
P4-P5 + 



[3|2|4> 



[3|5|4>. 



(4.4) 



(4.5) 



(4.6) 
(4.7) 



The result (J4.4J) is much more compact than the expression obtained from a Feynman diagram 
calculation, with which it agrees. See Section 5 for the Feynman results for this process in 
terms of massive spinor products. 



4.1 Results for Helicity Conserving Amplitudes 

Here we give all the helicity conserving QCD amplitudes for qq — > ggg. By helicity conserving 
we mean those amplitudes where the spin polarizations of the fermions are H — , in the sense 
described in Section 2. Whether these labels actually correspond physically to helicity depends 
on the choice of k . We choose to mark adjacent gluons, so that each amplitude has contributing 
recursive diagrams of the form of Fig. 2a and Fig. 2b, that is, we have a diagram with an 
internal fermion and a diagram with an internal gluon. The vanishing of the 2 — > 2 amplitude 
M + h+ simplifies those cases where there is a majority of positive helicity gluons. In particular, 
the diagrams with an internal gluon vanish. In the remaining cases, we evaluate such diagrams 
in the same way as in Ref. jHj, using identities such as 

IAP] - I*™, (4.8) 

(PB) = (4.9) 
\JP] 

with i and j as in ([3.1jh These identities hold only when the marked particles i and j are 
massless. 

The results presented here are valid for arbitrary spin polarizations. Choosing a polarization 
basis amounts to choosing the vector ko, and when this is done the expressions below will 
simplify. In the helicity basis for example, in which we choose ko to be parallel to the line of 
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approach of the fermions, the building block M + vanishes. This causes the first term in 
each of the mostly-minus amplitudes below to vanish. 



M+-+-+ = [5|1|4> 
where i — 3, j — 4 and 



m(14)(42)[53] + [15] (42)[3|2|4) - (14) [32] [5| 1 14) 

8 P5 ■ Pl P2 ■ P3 Pa ■ P5 (43) 

P2 - Pa 



(4.10) 



Pa ■ P5 = Pa ■ Pb 



[3|2|4> 



[3|5|4>. 



M + - ++ - = [4|l|5)x 
m[14][32](54) + [14] (42) [3 1 2 1 5) + 2p 5 ■ pi(15)[32] + m(15)(42)[43]" 



where i — 3, j — 4 and 



P4 -P5 



-P3 PA -P5 P5 'Pi (43) 



^•P5 + ^g[3|5|4), |4] = |4]- 



-2 p 2 • gg) 
[3|2|4) 



|3]. 



(4.11) 



M+~ ++ = [4|2|3)x 
m[42][15](43) - [42] ( 14) [5 1 1 1 3> - 2p 2 • p 3 [15](32) + m(14)(32)[54]~ 



&P2-P3 P3 - Pa P5 - Pi (54) 



where « = 5, j = 4 and 



P3 -Pa 



P3 ' P4+ [5W) [ 1 1 )? 



141 Ml (~ 2 Pi -P5) | K1 

|4] = |4] " [5|1|4) |5] - 



(4.12) 



m[21](4, c 



C\3 



[3|2|4) 



(34) [(35)2p 5 • Pl + (34) [4|2|5)j ( Pl + p 2 f 
m(42)(15)[43] - (42)[14][3|1|5) - 2p 2 • p 3 (15)[32] + m[14][32](54)' 



8p 2 -P3 P3 - Pa Pb -Pi[54] 



where i = 4, j = 5, and 



P3 -P4 = P3 -Pa + —\ i 4 l 3 l 5 )» 



[4|1|5) 



(42) = (42) 



(2gi • gg) 
[4|1|5) 



(52). 



(4.13) 



12 



m[21](53) z 



V2 (43) (45) [(53)2 p 2 ■ p 3 + (54) [4|2|3>] ( Pl + p 2 f 

m[14][42](53) + ( 15) [42] [4| 2 1 3> - [14] (32) [4j lj 5> " 



(4.14) 



+ [4|1|5> 



&P5-P1 P2-P3 P4 -P5[34] 



where i — 4, j — 3, and 



• /;.-, = /'i • /'.-, • 'J.J^ i 5 3). 



[4|2|3> 



M + " 



[5|1|4> x 



m[21] (43) J 



+ 



v/2 (45) [(53)2 P2-P3 + <54>[4|2|3>] ( Pl + p 2 y 
m(14)(32)[54] + (14)[42][5|2|3) + 2p 5 • pi[15](32) + m[15][42](43)' 



(4.15) 



&P2-P3P4- Pb Pb ■ Pi [34] 



where i — 4, j — 3, and 



. P2-P3 

P4-p 5 = P4-P5+ [4| 2 |3) ^I 5 I 3 /' 

14) = |4) + 2P2 ' P3 I3) 
1 ; 1 ; [4|2|3) 1 ; 



_ = m (54)[4|2|3)[12][45] 

4 [45]> 5 .p lP! -p 3 [34] 1 ' 

where 

< 5? > = < 54 > + #jf (53> - 

M+-+++ = (4.17) 

The amplitudes with fermion helicities — h can be obtained from those above by complex 
conjugation. 

4.2 Results for Helicity Flip Amplitudes 

We now consider the helicity flip amplitudes. These have fermion spin polarization labels ±±. 
In this case we are unable to evaluate diagrams with an internal gluon due to products such as 
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(Pk) where k is massive. For those amplitudes in which all gluons have the same helicity, the 
internal gluon diagrams vanish anyway: 



M++ = 

where i — 4, j — 5 and 



m(43) 



m(15)(42)[43] + [14] (32) [4| 1 1 5) - [14] (42) [3 1 1 1 5> " 



-4 p 5 - Pl p 2 -p 3 [34] 2 [54] 



(4.18) 



|4) = |4> + 



2 p 5 ■ Pi 
[4|1|5> 



|5). 



M +++++ = m[54] 
where i = 3, j = 4 and 



m(14)(32)(54) + (14) (42) [3 1 2 1 5) - (15)(42)[3|2|4) 
-4p 5 - Pl p 2 -p 3 {45) 2 {43) 

|41 = |41 + 2P2 - P3 |31 
11 11 [3|2|4> 1 J 



(4.19) 



The amplitudes M +++ and M are obtained from those above by complex conjugation. 

For the remaining amplitudes we resort to Feynman diagrams. 

5 Feynman Results 

Here we give results for qq — > ggg derived from Feynman rules. Note that in a given amplitude 
all the helicities can be flipped by complex conjugation. 



M + ~ 



-+- 



= -[4|2|3> 
+ (35) 



m[14][42](53) + ( 15) [42] [4| 2 1 3) - [14] (32) [4 1 1 1 5) 
8 p 5 ■ Pi P2 ■ P3 Ps ■ Pa [54] 

m[14][42](53) + (15) [42] [4| 2 1 3> - [14] (32) [4j 1 1 5> 



+ 



(35> s 



Sp 2 -P3P3- Pa Pa • P5 



[14] (32) + (13) [42] [14] (52) + (15) [42] 



(34)(45)( Pl +p 2 ) 



[54] 



+ 



[34] 



(5.1) 



-++- 



-[4|1|5> 
+ [43] [4|1|5> 
- [43] 



-m(15)(52)[43] + (15) [32] [4| 1 15> - [14] (52) [3 1 2 1 5) 
8p5-PiP2- P3 Pa ■ P5 (53) 

[14](52) + (15)[42] 



4 p b - pi p 3 -p 4 (53)[54]_ 
[14](52)[3]lj5> + (15)[32][4|1|5) - m(15)(52)[43] 



[43] 2 (35) 



4 Pb ■ Pi P3 ■ Pa Pa • P5 



[14] (52) + (15) [42] [13] (52) + (15) [32] 



2 p 3 -p 4 [54](pi +p 2 f 



(53) 



+ 



(54)] 



(5.2) 



14 



= [3|2|4> 
+ (45)[3|2|4) 
+ (45) 



m[13][32](54) - [13] (42) [3 1 1 15> + (15) [32] [3| 2 14> 



-8 p 5 ■ pi p 2 ■ P3 p 3 ■ Pa [53] 
[13] (42) + (14) [32] 



4 p 2 -p 3 p A -p 5 (43)[35] 
[13] (42) [3|2|5) + (15) [32] [3|2|4) + m[13] [32] (54)] 



+ 



(45) 2 [53] 



8 P2 ■ P3 P3 ■ P4 Pi ■ Pb 

■ [13] (52) + (15) [32] [13] (42) + (14) [32] 



2p 4 -p 5 (43}( Pl +p 2 y 



[43] 



[53] 



(5.3) 



The corresponding helicity flip amplitudes can be obtained from these simply by altering 
the types of brackets. For example, suppose we wish to extract M 1 from M H h ~ given 
above. We can achieve this by changing brackets as follows: 



[lk] 
(lk) 



(lk), 
(lk), 



(5.4) 
(5.5) 



where k is massless. Sandwich products such as [4|1|5) are not changed. This transformation 
results in 



M~ 



= -[4|2|3) 
+ (35) 



m(14)[42](53) + (15) [42] [4|2|3) - ( 14) (32) [4 1 1 1 5> 
8 p 5 ■ Pi P2 ■ Pz P3 ■ Pi [54] 

m(14)[42](53) + (15) [42] [4|2|3> - ( 14) (32) [4| 1 1 5> 



+ 



(35) 1 



8 P2 ■ P3 P3 ■ P4 Pi ■ Pb 

(14) (32) + (13) [42] (14) (52) + (15) [42] 



(34)(45)( Pl +p 2 ) 



[54] 



[34] 



(5.6) 



Other amplitudes can be found by analogous bracket alterations. 



6 Summary 

We have calculated all the partial spin amplitudes for the qq — > ggg scattering process where q is 
a massive fermion. For most of the partial amplitudes we were able to use the BCFW recursion 
relations to obtain fairly compact expressions. This was achieved by following the idea of 
Ref. of stripping lower point amplitudes of their external fermion wavefunctions before 
inserting them into the recursion. We used a particular representation of massive spinors, 
along the lines of the appendix of Ref. [2E], to define massive spinor products. In this method 
information regarding the polarization of the fermion spins is contained in the definition of the 
spinor products, rather than explicitly in the amplitude. 
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We derived very compact results for the helicity conserving partial amplitudes. This sim- 
plicity can be attributed to the vanishing of certain 2 — > 2 scattering amplitudes, which meant 
that in some cases we were able to avoid diagrams with internal gluons. We were unable to 
treat the helicity flip amplitudes in the same way (except for the case where all the gluon helic- 
ities are the same), since we were unable to evaluate the corresponding recursive diagrams with 
internal gluons. For these amplitudes we instead provided expressions derived from Feynman 
diagrams, also in terms of massive spinor products. 

These results represent an interesting test of the BCFW recursion relations HDJ. The 
massive spinor products we used are well suited to calculations with massive fermions. Appli- 
cation of these techniques to higher order processes with massive fermions, such as qq — > gggg, 
should be possible though would be accompanied by an increase in complexity. 
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Appendices 



A Notation and Conventions 

We have used products of Dirac spinors 

[ij]=u + (i)u-(j), (ij) =u~(i)u + {]) ) (i ] )=u ± {i)u ± {j) (A.l) 
with massive Pi,Pj. To evaluate these we use two arbitrary four vectors k and hi, such that 

kl = 0, k\ = -1, k • h = 0. (A.2) 

Then 

r . (Pi • h) {pj ■ h) - {pj ■ ko) {pi ■ h) - ie^vpcKpXp] K / k 

VlPi • k o){Pj ■ ko) 
(Pj ■ h)(Pi ■ h) ~ (Pj ■ h)(pi ■ k ) - ievvpcKp^fil 

VlPi ' k o){Pj ■ ko) 

(ij) = mi( P3 f° ) (A.5) 
\Pi ■ ko) 

where mj is negative if i is an antiparticle. Different choices of ko correspond to different choices 
of the quantization axis of a massive fermion's spin, as described in Section 2. 

We use the notation 

u + (i) f u + (j) = \i\p\j) = [ip](pj) + (ip)(pj), (A.6) 
u~(i)fu~(j) = (i\p\j] = (ip}[pj} + (ip)(pj). (A.7) 

Whereas for massless vectors k i: kj we have the familiar relation 2k, r kj = in the massive 
case this is extended to 

2p i 'P j = (ij)\ji] + (ij) 2 . (A.8) 
For any massive i,j and massless k, I we have 

(ik)(jl) = (il)(jk), (A.9) 
(ik)[li] + [ik](li) = m^lk], (A. 10) 

(Pfdj/iv^to) = (A.ll) 



u 



± 



The Schouten identity holds, 

(a b)(c d) + (a c)(d b) + (a d)(b c) = 0. (A.12) 
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For gluon polarization vectors we use 



4m = s ' ( *>y; (p) . (A.i3) 

V2 {kp) 

e -( P ,k) = (A.14) 



which take the slashed form 



= V2 M+( ^ +(p) ( + ) M - (p) "- (fc) ) (A.15) 

re?,*) = ^" + (p)"4fc)+^-^-(p). (A.i6) 

We use a shorthand form for the amplitude in which we display only the helicities of the particles 
involved. So for example, 

M(q+,qz ; 3 + ,4~,5 + ) ~ M {+ - + - +) . (A.17) 



B Colour Decomposition 

The calculation of multi-parton scattering amplitudes in perturbative QCD becomes problem- 
atic very quickly as the number of partons increases, due to the sheer number of diagrams and 
the complicated gluon self-interactions. One technique to circumvent this is to split the set 
of all Feynman diagrams contributing to a particular amplitude into gauge invariant subsets. 
Then different gauges can be used in the evaluation of each subset. This simplifies the overall 
calculation considerably. Each subset of Feynman diagrams is called a partial amplitude. We 
use the normalization 

Ty(T a T b ) = 5 a b- (B.l) 

where the T l are matrices of the fundamental representation of £77(3). This convention leads 
to colour-ordered Feynman rules as given in [31] ■ The colour decomposition for processes with 
a pair of quarks is then 

A(q,q; 9l ,g 2 ■ ■ ■ g n ) = ^ • • • T<Mn) ) ^{q, q ; <?,(!), g*{2) • • • 9*(n)). (B.2) 

Here a is the set of all distinct cyclic orderings of the gluons. The colour information in a given 
amplitude is contained purely in the group theoretical prefactors, while all the kinematical 
information is contained in the partial amplitudes M(a). It is useful to note that amplitudes 
in QED can be obtained from amplitudes in QCD by replacing all the colour matrices T A with 
the identity matrix. For references and a more detailed description of colour decomposition, 
the reader is directed to Ref. |5T] . 
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